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Two types of geometrical asymmetries have been analyzed numeri- 
cally to determine their influence on the externally visible caustic 
patterns from a drawdown zone. The first of these asymmetries involved 
axial displacement of identical opposing profiles of the drawdown zone, 
while the second involved opposing profiles of differing slope positioned 
with their inflection points in axial alignment. The results of these 
studies are presented graphically as an aid to interpreting the asym- 
metries of actual caustic patterns generated by illuminating real 
drawdown zone samples in optical fiber drawing. 

I. INTRODUCTION 

In the present study, we demonstrate how deviations from rotational 
symmetry in the melt zone give rise to asymmetries in the resulting 
caustic pattern. In Ref. 1, rotationally symmetric profiles of melt zones 
were generated by averaging measurements of data taken from solidified 
melt zone samples. These data provided symmetrically idealized data 
sets for the study of the two principal types of caustics as functions of 
the rate of change in the cross section. Changes in the maximum rate of 
change in the cross section were simulated by scaling* the data sets, and 
the analysis was carried out using the algorithms developed earlier in 
Ref. 2. These algorithms describe families of 2- and 3-intercept caustics 
as shown in Fig. la, where the scale factor is unity. The incident angle, 
/3, for an axial light ray at the inflection point, /, is 52.9 degrees on both 
sides of the cross section. However, in Ref. 2 it was shown that the 
samples were usually asymmetric, resulting in pronounced asymmetric 
caustic patterns, especially when the caustic rays emerged nearly tangent 
to the surface. 



* That is, by scaling the axial and transverse data values differently with the scale factor, 
SF = (axial scale)/(transverse scale). 
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Fig. 1 — Drawdown profiles for a melt zone with an average 7 of 1.323 and a preform 
radius, R, of 0.1217 in. (3.09 mm). The data scale factor, SF, is 1. The rays labeled D rep- 
resent the 3-intercept caustic trajectories and the rays labeled 1 represent the first 
2-intercept caustic trajectories for the following cases: (a) Symmetric drawdown profiles 
with 7 = 1.323 (0 = 52.9 degrees) on both sides and no shift (A/yR = 0.0, by = percent), 
(b) Asymmetric drawdown profiles with 7 = 1.323 (/3 = 52.9 degrees) on both sides and 
an axial shift of A = 0.05-in. (1.27 mm) of the left side ahead of the right side (A/yR = 0.311, 
AIR = ±0.41). Here by = percent, (c) Asymmetric drawdown profiles with 7 = 1.455 (/3 
= 55.5 degrees) on the left side and 7 = 1.1907 (/3 = 50.0 degrees) on the right side (average 
7 = 1.323, by = ±10 percent). Here A/yR = 0.0. 

Two types of asymmetric drawdown zones were studied, based upon 
the symmetric data set generated from the fourth sample. 1,2 In the 
simplest case, the asymmetry was introduced by axially shifting the 
profile data for one side of the cross section with respect to the other. An 
asymmetric cross section with the caustic rays is shown in Fig. lb, where 
the relative axial shift between the two profiles is 41 percent of the 
preform radius* and the scale factor is 1.0. The asymmetry produces 
significant counterclockwise rotations in the propagation directions of 
the emerging caustic rays when viewed with the downstream profile as 
shown on the left. 

A second type of asymmetry was produced by scaling the opposite 
sides differently. In this study, the inflection points of the opposite 
(differently scaled) profiles were kept in axial alignment so that the re- 
sulting shift in the caustic ray pattern would be due solely to the dif- 
ference in the slopes, as shown in Fig. lc. There the scale factors of op- 
posite sides are 0.9 and 1.1. These produce counterclockwise rotations 
of the emerging caustic rays comparable to those shown in Fig. lb and 
toward the elongated profile as plotted. 

In addition, for both types of asymmetry, the overall scaling was 
changed so that the effects of these asymmetries could be assessed for 
greater or lesser average gradients than those of the original data. These 
scale factors and associated average slopes, 7, and incident angle values, 
Pa, are tabulated in Table I. Here 7 is the average of the slopes of the two 
outer normals at the inflection points, dx/dyj, and Pa = —arc tan 7. 

The next two sections present the results from the computer simula- 
tions including diagrams of significant cases for both types of asymme- 

* R =0.1217 in. or 3.09 mm. 
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Table I — Scale factors and average slopes and angles 





Average Outer 


Average Incidence 


Scale Factor 


Normal Slope y 


Angle 


A (degrees) 


1.300 


1.720 




59.8 


1.200 


1.588 




57.8 


1.000 


1.323 




52.9 


0.775 


1.025 




45.7 


0.700 


0.926 




42.8 


0.520 


0.688 




34.5 


0.450 


0.595 




30.8 


0.340 


0.450 




24.2 



tries. The reader who requires a less detailed synopsis of these results 
is referred to the discussion portions of these sections and to Section V, 
Summary. 

II. SHIFT ASYMMETRY 

The effects of axial misalignments in the profiles of the drawdown zone 
were studied by computing the ray paths for examples with the profiles 
of opposite sides shifted various amounts, A. The results of these com- 
putations provide a complete description of the caustic patterns for an 
asymmetric drawdown zone with a range of profile shifts for each y given 
in Table I. 

To compare the varied scalings, an asymmetry parameter, A/yR, was 
established, where R is the preform radius. The asymmetry was varied 
from < | A/yR \ < 1.20, which represents a broader range than is ever 
likely to actually occur in fiber drawing. The specific incremental changes 
in asymmetry were chosen to illustrate what occurs as the asymmetry 
increases. 

An extreme example of the effects of this type of asymmetry can be 
seen by comparing Fig. 2a with Fig. 2b, which has the same average slope 
but a large asymmetry. Emergent 2-intercept caustics are labeled nu- 
merically to indicate which of the family they represent; the 3-intercept 
caustics are labeled D. The unlabeled rays represent caustic rays of the 
1 -intercept family which emerge at the inflection point for 7s less than 
0.939. In Fig. 2b, the 3-intercept caustics have rotated counterclockwise, 
or in a positive direction, toward the leading side. Extra 2-intercept 
caustics (labeled 2', 2", 3') have merged on the leading side, while all 
2-intercept caustics have disappeared from the trailing side. Those 
2-intercept caustics on the +A side have rotated clockwise, in a negative 
direction, away from the leading side. 

2.1 Graphical results 

Figures 3a to 3d are simplified versions of the format used in Figs. 1 
and 2 which show both the internal ray paths and the external caustic 
trajectories. In these simplified diagrams, the internal ray paths have 
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Fig. 2 — Drawdown profiles for a melt zone with an average 7 of 0.688 (SF = 0.52). The 
rays labeled D represent the 3- intercept caustic trajectories, while the rays labeled 1, 2, 
3, etc., represent the first, second, third, etc., 2-intercept caustic trajectories for the fol- 
lowing cases; (a) Symmetric drawdown profiles with 7 = 0.688 (/3 = 34.5 degrees) on both 
sides and no shift, A/yR = 0.0, f>y = percent, (b) Asymmetric drawdown profiles with 
7 = 0.688 (0 = 34.5 degrees) on both sides and a shift, A = 0.075 in. (1.905 mm) of the left 
side ahead of the right side {A/yR = ±0.896, A/R = ±0.616). Here, 6y = percent, 
(c) Asymmetric drawdown profiles with 7 = 0.894 (/3 = 41.8 degrees) on the left side and 
7 = 0.482 (0 = 25.7 degrees) on the right side (7 = 0.688, 67 = ±30 percent). Here, 
A/yR = 0.0. 

been eliminated, and the drawdown profiles have been reduced to tan- 
gent lines at the inflection points. The profiles, shown as thick lines, meet 
in the symmetric case and progressively move farther apart to depict 
increasing ±As. The external caustic trajectories are shown schematically 
as rays radiating from the point A = 0. The 3-intercept caustics are again 
labeled D and, in the unusual case of a second 3-intercept caustic, D2 
is used. Whenever the ray which internally reflects from the inflection 
point, /, is refracted from the opposite side without forming an external 
caustic, it is shown by an unlabeled dashed line. In Figs. 3a to 3d, the 
diagrams are drawn with the leading side as the left profile, consistent 
with Figs, lb and 2b. We then define a positive, or forward, rotation of 
caustic angles to be counterclockwise, or toward the leading side. 

We begin with gradually tapered profiles that emit little light. As the 
taper is increased, we observe an increase in the amount of light emitted 
and in the complexity of the caustic structure. The first case, with 7 = 
1.72, is shown in Fig. 3a. Both caustics present initially rotate in a forward 
direction as the asymmetry increases. However, the 3-intercept caustics 
quickly disappear from the trailing side by internal reflection. The 2- 
intercept caustics continue to rotate in a positive direction on the trailing 
side, but on the leading side they reverse direction and rotate in a neg- 
ative direction at an asymmetry A/yR of about 0.40. The situation is 
similar for 7 = 1.588, not shown. The 3-intercept caustics also rotate 
forward. Since the emergent angle for the symmetric case is smaller this 
time (163.48 degrees as opposed to 172.66 degrees), the caustic does not 
internally reflect as soon. The 2-intercept caustic initially rotates for- 
ward. This time, the direction reversal on the leading side occurs earlier 
at a A/yR of 0.33 and is more pronounced. Reducing 7 to 1.323 does not 
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Fig. 3 — Diagrams of external caustic trajectories for drawdowns of varying asymmetry 
due to axial shifts, A/7/?, of varying amounts as follows: (a) For 7 = 1.720 (0 = 59.8 degrees), 
A/7/? = ±0.000, ±0.118, ±0.236, ±0.354, ±0.472, ±0.596. (b) For 7 = 1.025, (0 = 45.7 de- 
grees), A/7/? = ±0.000, ±0.199, ±0.410, ±0.596, ±0.801. (c) For 7 = 0.926 {0 = 42.8 degrees); 
A/7/? = ±0.000, ±0.223, ±0.441, ±0.664, ±0.888. (d) For 7 = 0.450 (0 = 24.2 degrees); 
A/7/? = ±0.000, ±0.183, ±0.365, ±0.548, ±0.730, ±0.913. Here the sides of the melt zone 
are represented by heavy lines sloped (1 degrees below horizontal and shifted vertically 
in proportion to the axial shift, A. The external caustic trajectories are shown as rays em- 
anating from a common center at their appropriate angles and labeled with D or D2 to 
indicate a first or second 3-intercept caustic and 1, 2, or 3, etc. to indicate a first, second 
or third, etc. 2-intercept caustic. Dashed lines represent emergent / rays which do not form 
caustics. 
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alter the situation significantly, except that both caustics are present 
for the entire range of asymmetry. In all cases, the rotation of all caustics 
on the trailing side has been consistently positive. On the leading side, 
the 2-intercept caustic's forward rotation reverses at a A/7/? of only 0.25. 
At the same time, a reversal in the forward rotation of the 3-intercept 
caustic also becomes apparent at an asymmetry of 0.75. 

Figure 3b shows the case of 7 = 1.025, where (3 = 45.7 degrees, i.e., very 
close to the 2-intercept caustic's extinction angle. This caustic emerges 
at a minimum angle, 0, for the symmetric case, and progressively forms 
larger angles as A increases; i.e., it rotates negatively on the leading side 
and positively on the trailing side. As discussed in Part I, 1 the rays which 
are only incident once on the first side form a catacaustic at the inflection 
point. If the catacaustic ray then refracts from the second side, it is 
usually seen as the first 2-intercept caustic ray. However, the dashed 
lines shown in the diagrams for the smaller asymmetries in Fig. 3b depict 
emergent rays which form internal catacaustics at the inflection points 
but are not seen as external caustics. We call these rays the / rays. Their 
angles are at least half a degree different from the 2-intercept caustic 
angles. Once the asymmetry has increased to a A/yR of 0.80, however, 
these two rays have recombined. The 3-intercept caustic rotates in the 
forward direction over the asymmetry range studied. 

Continuing further, the symmetric case with 7 = 0.926 has no emer- 
gent 2-intercept caustic (see Fig. 3c). However, once the asymmetry has 
increased to A/yR = 0.441, the original 2-intercept caustic has re- 
emerged and subsequently rotates to larger 6 angles on both sides. Again, 
the (noncaustic) / ray is initially distinct from the 2-intercept caustic 
ray, and they merge as the asymmetry increases. It should be noted that 
the / ray initially emerges on the forward side without forming a caustic 
before the 2-intercept caustic emerges. The 3-intercept caustic continues 
to rotate in a forward direction up to the largest calculated asymmetry, 
at which point it reverses direction on the forward side. 

When 7 is 0.688, there are three 2-intercept caustics for the symmetric 
case caused by multiple-folding in the fan of emerging light. This is the 
case shown earlier in Fig. 2a (A/yR = 0.0) and Fig. 2b {A/yR = 0.896). 
The folding geometry is illustrated in the left and center diagrams in 
Fig. 4a. Here the caustic rays are labeled 1, 2, 3, etc., consistent with the 
present convention, while the bounding rays are labeled LI, L2, L3, etc. 
The lower L numbers represent bounding rays originally propagating 
nearest the surface of the sample, with successive higher numbers pro- 
ceeding inward. The left diagram in Fig. 4a shows that the fan of light 
which lies between bounding rays L3 and L4 is folded to form the three 
2-intercept caustics with the first and third very nearly superimposed. 
As the asymmetry is increased, most of the 2-intercept caustics disappear 
on the trailing side while several additional ones appear on the leading 
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side. All these caustics rotate in a negative direction except for the 
original 2-intercept caustic on the leading side, which remains stationary. 
This unusual circumstance occurs because the / ray does not form an 
external caustic until its point of emergence has moved upstream of the 
melt zone to where the radius is constant. 

In Fig. 2b, we can see that, since the profile on the initial incident side 
does not change, the internal trajectory of the / ray remains the same. 
Consequently, as long as caustic ray 1 refracts from the leading side its 
external angle will be constant, regardless of the magnitude of A. On the 
trailing side, the internal / ray is incident along that portion of the melt 
zone profile for which the radius is changing very rapidly; therefore, al- 
though it emerges for asymmetries of less than 0.3 it does not form an 
external caustic. However, this / ray reappears along with a first 
2-intercept caustic when the asymmetry A/yR reaches 1.2. In addition, 
a new 2-intercept caustic emerges on the leading side. This fourth 2- 
intercept caustic is associated with the closing of the break between the 
L1-L2 fan of light and the L3-L4 fan of light, illustrated in the center 
diagram in Fig. 4a. Consequently, this new caustic is folded downstream 
in the same way as the original first 2-intercept caustic. It originates from 
a coaxial ray initially propagating down the sample at a greater radius 
than the second 2-intercept caustic, but well inside the ray which origi- 
nates the 3-intercept caustic, as shown in Fig. 2b. 

Several other minor 2-intercept caustics appear and are labeled 2', 2", 
and 3' in Fig. 4a center. These lie very close — in terms of both originating 
axial rays and emergent angles — to the second or third 2-intercept 
caustics. Each is associated with another fold in the fan of light.* The 
3-intercept caustics in this case (7 = 0.688) initially rotate positively. 
Eventually, on the leading side this caustic reverses its rotation at 
A/yR ~ 0.6 and then extinguishes at an asymmetry of about 0.8. On the 
trailing side, it continues to rotate in a forward direction. The additional 
2-intercept caustics also appear in the case where 7 = 0.595. Again the 
first 2-intercept caustic, when it emerges, remains stationary. The others 
in this family all rotate in the negative direction. The first, second, and 
third 2-intercept caustics emerge on the leading side for all asymmetries 
and are extinguished when the asymmetry is between 0.83 and 1.1 on 
the trailing side. Minor caustics associated with the second and third 
major caustics were found on the trailing side, but not on the leading side. 
The fourth caustic emerges at greater asymmetry than in the previous 
case of y = 0.688 and again only on the leading side. This case, with 7 
= 0.595, is of particular interest because it exhibits for the first time a 
second 3-intercept caustic. The new caustic is associated with a folding 



* It should be noted that, since the algorithm computes the ray trajectories for an equally 
spaced sequence of axial rays, finding these minor perturbations in emergent ray angles 
depends somewhat on the ray spacing chosen. 
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Fig. 4 — (a) Schematic diagrams showing the external limiting rays for the 2-intercept 
light paths. The drawdown melt zone profiles are omitted, but the orientation is the same 
as in Figs. 1 and 2. The vertical arrow represents the incoming beam of collimated light 
illuminating the sample which in turn deflects the light into the plotted ray paths in the 
far field. Here, the solid lines represent the caustic trajectories and the broken lines rep- 
resent the bounding rays limited by internal reflection (see Ref. 1). The circumferential 
arrows follow the continuous fan of light from the outermost bounding ray to the innermost 
bounding ray including all caustic rays. The bounding rays are numbered in sequence, 
beginning with LI as the ray which originates nearest the surface of the sample and in- 
creasing inward L2, L3, etc. The caustic rays themselves are numbered 1, 2, 3, etc. to in- 
dicate that they represent the first, second, or third, etc. of the 2-intercept caustics. The 
asymmetries are as follows: Right — Symmetric, with 7 = 0.688 (0 = 34.5 degrees) on both 
sides and A/yR = 0.0, by = percent. Center — Asymmetric with 7 = 0.688 (0 = 34.5 de- 
grees) on both sides and a shift A = 0.075 in. (1.905 mm) of the leftside ahead of the right 
side (A/yR = ±0.896, MR = ±0.616). Here, 67 = percent. Left— Asymmetric with 7 = 
0.894 {0 = 41.8 degrees) on the left side and 7 = 0.482 (0 = 25.7 degrees) on the right side 
(7 = 0.688, 67 = ±30 percent). Here, A = 0.0. (b) Schematic diagram showing external 
limiting rays for the 3-intercept light paths. The same ray identifications apply to the di- 
agram as apply to the diagrams in Fig. 4a except that the caustics are numbered D and 
D2, representing the original and second 3-intercept caustics. The asymmetry is given 
by 7 = 0.595 (0 = 30.8 degrees) on both sides and a shift, A = 0.08 in. (2.032 mm) of the 
left side ahead of the right side ( A/yR = ±1.105, A/R = ±0.657). Here, 67 = percent. 

in the fan of light in the direction opposite that of the original 3-intercept 
caustic. It is formed from a ray initially propagating down the fiber at 
a radius smaller than that ray which forms the first 3-intercept caustic. 
The folding pattern is shown in Fig. 4b, where the new ray is labeled D2. 
The new caustic emerges only for large asymmetry and only on the 
leading side. Meanwhile, the original 3-intercept caustic, as before, ro- 
tates initially in a positive direction but eventually reverses with in- 
creasing asymmetry. 

The last case studied, 7 = 0.450, is shown in Fig. 3d. Here the first 
three 2-intercept caustics emerge over the whole range of asymmetries. 
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The second and third caustics rotate in a negative direction; the first is 
again stationary. The fourth caustic emerges on the leading side and only 
for the largest asymmetry. The first 3-intercept caustic appears only on 
the trailing side and rotates downstream with increasing asymmetry. 
The second 3-intercept caustic appears alone on the leading side and only 
for larger asymmetries. 

2.2 Discussion 

The diagrams in Fig. 3 illustrate the far-field caustic response to in- 
creasing axial asymmetry A/yR for fixed 7. We now present a general 
quantitative description of the caustic behavior due to changing geom- 
etry. Figure 5 shows the propagation angles, 0, of the major caustics as 
functions of A/yR for various 7. The values, 0_, for the trailing side are 
plotted to the left and the values, 0+, for the leading side are plotted to 
the right of the zero asymmetry line. Figures 5a and 5b show the first 
2-intercept and the 3-intercept caustics, respectively. The behavior of 
the remaining major 2-intercept caustics, the second, third, and fourth, 
are plotted in Fig. 5c. Each curve on these plots represents a fixed outer 
normal slope y over the range of asymmetry tested. Those curves with 
arrowheads at the ends indicate caustics which emerge for higher 
asymmetry but were not calculated. The absence of an arrowhead in- 
dicates that the caustic ceases at or just beyond that point. 

Other measures of caustic behavior are given in Fig. 6, for cases where 
the caustics emerge on both sides of the drawdown. In terms of actual 
three-dimensional drawdown samples, the caustic trajectories form cones 
of light whose total included angle is 0- + 0+. The average sum, 0s = l k(Q- 
+ 6+), is plotted to show the half cone angle as a function of asymmetry. 
The average difference, do = %(0- — 0+), is plotted to show the rotation 
of the cone as a function of asymmetry. Figures 6a and 6c refer to various 
2-intercept caustics, while Fig. 6b refers to the 3-intercept caustics. 

Consider the first 2-intercept caustic shown in Figs. 5a and 6a. In Fig. 
5a, starting with large 7, on the leading side 6 U1 decreases, reaches a 
minimum, and slightly increases as A/yR increases. As the profile be- 
comes blunter, the minimum becomes more pronounced and moves 
toward the symmetric case, A/yR -* 0. When 7 = 0.926, the drop in Q m 
as A/7/? — becomes so large that the caustic is extinguished for small 
asymmetries and only the ends of the curve remain. This means that, 
although the first 2-intercept caustic is extinguished by internal re- 
flection in the symmetric sample of 7 = 0.926, the development of suf- 
ficient asymmetry causes an almost simultaneous reappearance on both 
sides of the sample. Figure 6a also shows that the average angle, S U } of 
the first 2-intercept caustic decreases as A/yR — 0. This decrease is 
greatest for profiles with the lowest 7s. In other words, the caustic cone 
angle increases as the asymmetry increases, especially in blunt samples. 
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Fig. 5 — Plots of the 6 angles of the caustic trajectories for various 7s as functions of the 
axial shift, A. Here the values of the angles, 6+, of the caustics emergent on the leading side 
are plotted on the right and the values of the caustic angles, 6-, of those emergent on the 
trailing side are plotted on the left. These results are arranged as follows: (a) The first 
2-intercept caustic angles, designated 8 UI , and the / ray propagation angle, 6 1 , vs A/7/?, 
(b) The first and second 3-intercept caustic angles, 0° and 0" 2 , vs A/7/?, (c) The second, 
third, and fourth 2-intercept caustic angles, d" 2 , U3 and UA , vs A/7/?. 

At the same time, the dp 1 vs A/yR curves in Fig. 6a show that the first 
2-intercept caustic cone, which is initially symmetric about the drawing 
axis, first rotates toward the leading side, reverses, and then rotates back 
toward the axis as the asymmetry increases. 

The 3-intercept caustic exhibits similar behavior, as shown in Figs. 
5b and 6b, although there are significant differences. Again, there is a 
minimum in the caustic angle vs asymmetry curve that appears on the 
leading profile side of the 6 D vs A/yR plot. However, the minimum does 
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Fig. 6— Plots of the averaged sums, Q s = V2C- + "+). and the averaged differences, fl D 
= y 2 (fl_ - 0+), vs the magnitude of the axial shift | A\/yR for various 7s. These results are 
arranged as shown in 6a, 6b, and 6c. (a) The first 2-intercept caustic angles, Q m , averaged 
sums and differences vs | l\/yR. (Figs. 6b and 6c on following pages.) 

not move toward the symmetry axis as 7 decreases. Also, except in the 
region near the extinction points, the D vs AAyfl plots are morphologic- 
ally similar. In Fig. 5b, the magnitude of the depression in the curve still 
grows as 7 decreases until, by 7 = 0.450, the caustic is extinguished on 
the leading side. Figure 6b shows that the caustic cone half-angles, 8$, 
are less influenced by the asymmetry than they were for the 2-intercept 
caustics (Fig. 6a). In addition, the 3-intercept caustic cone rotates more 
rapidly toward the leading side before it reverses direction, as shown by 
the 6% curves in Fig. 6b. The behavior of the second, third, and fourth 
2-intercept caustics is illustrated in Figs. 5c and 6c. These caustics are 
distinct because their directions of rotation are opposite those of the first 
two major caustics. This backward rotation, with increasing A/7/?, of 
the second and third 2-intercept caustics parallels the reverse 6a rotation 
with increasing 7 noted in the earlier "symmetric" study. 1 At the same 
time, the 6s values in the second and third 2-intercept caustics remain 
almost constant, indicating that the shapes of these caustic cones are 
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Fig. 6(b) — The 3-intercept caustic angles, d D , averaged sums and differences 
vs|A|/ 7 fl. 



almost independent of the asymmetry, except near the extinction 
points. 

III. DIFFERENTIAL SLOPE ASYMMETRY 

The effects of asymmetric profiles of differing slopes were studied by 
changing the scaling of opposite sides of the profile by equal positive and 
negative percentages. These changes produce equal percentage changes 
in the slopes of the outer normals to the profiles, dx/dy. As shown in Part 
I, 1 the outer normal slope defines a unique incidence angle and is the 
most significant single parameter controlling the formation of the 
caustics. Therefore, we can define asymmetry by a ± change in 7, 
i.e., fry. 

An extreme example of the effects of this type of asymmetry is found 
in the comparison between Figs. 2a and 2c where in the latter dy = ±30 
percent. It can be seen that the external trajectories of the 3-intercept 
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Fig. 6(c) — The second and third 2-intercept caustic angles, U2 and 6 U '\ averaged sums 
and differences vs | A|/-y/?- 

caustics have rotated in the direction of increasing y, i.e., downstream 
on the shortened side and upstream on the lengthened side, which will 
be referred to as a forward (positive) rotation. At the same time, the 
external 2-intercept caustics have disappeared from the side of dimin- 
ished 7, while an additional 2-intercept caustic, the fourth, appears on 
the opposite side. Like the 3-intercept caustic rays, the first and third 
2-intercept caustic rays have rotated forward, while the direction of the 
second ray is unchanged. Indeed, only the insignificant 1 -intercept 
caustics exhibit a backward or negative rotation. A better understanding 
of this behavior can be obtained only in the context of a range of changing 
asymmetries for the various average 7s in Table I, which are discussed 
below. 

3.1 Graphical results 

A graphical description of the effects of this asymmetry covering a 
range of typical average 7s is presented in Fig. 7. As in Fig. 3, these dia- 
grams are simplified versions of the earlier formats showing only the 
external caustic trajectories. In most cases, the range of asymmetry, 8y, 
is from ±0 to ±30 percent. These represent a much greater range of 
asymmetry than might be expected to arise with actual samples. The 
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Fig. 7 — Diagrams of external caustic trajectories for drawdowns of varying asymmetry 
due to percentage differences in 7, 67, of varying amounts as follows: (a) 7 = 1.720 ± 0%, 
± 5%, ± 10%, ± 20%, ± 30%, ± 40%. (b) 7 = 1.025 ± 0%, ± 1%, ± 10%, ±20%, ± 30%. 
(c) 7 = 0.688 ± 0%, ± 10%, ± 20%, ± 25%, ± 30%, ± 35%. Here the sides of the melt zone 
are represented by heavy lines sloped degrees below horizontal. The external caustic 
trajectories are shown as rays emanating from a common center at their appropriate angles 
and labeled D to indicate a 3-intercept caustic and 1, 2, or 3, etc. to indicate a first, second, 
or third, etc. 2-intercept caustic. Dashed lines represent emergent / rays which do not form 
caustics. 
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specific incremental changes in asymmetry shown are not uniform, but 
were chosen to illustrate what occurs as the asymmetry increases. In 
every case, the diagram is drawn with the side of reduced 7 to the left 
so that a forward (positive) rotation of caustic rays is consistently 
counterclockwise on the diagram. In general, on the decreasing 7 side, 
increases and on the increasing 7 side, 6 decreases with increasing 
asymmetry. 

The first two cases, for which 7 = 1.720 (Fig. 7a) and 1.588, are similar; 
the rotation of all the caustic rays is forward and the 2-intercept caustic 
on the right-hand, or decreasing 7, side disappears by internal reflection. 
It initially propagates at a greater angle for 7 = 1.720 than it does for 
7 = 1.588, and its rotation due to increasing asymmetry is positive. 
Therefore, the right-hand 2-intercept caustic extinguishes at a 8y of only 
±10 percent for the case of 7 = 1.720, but must reach a 67 of ±14 percent 
before it extinguishes for 7 = 1.588. Decreasing 7 to 1.323 reduces 6 for 
all trajectories and permits forward rotations of both caustic rays on both 
sides for asymmetries up to 67 = ±30 percent without loss by internal 
reflection. 

Internal reflection can also cause a caustic ray to disappear when its 
propagation direction swings too far upstream (0 decreasing). This is 
shown in Fig. 7b where the 2-intercept caustic ray on the increased 7 side 
disappears as soon as the asymmetry begins to develop, because it 
emerges near the extinction region. 1 Very little asymmetry is then re- 
quired to rotate the left-hand ray below the critical angle. 

The case of 7 = 0.926 lies within the range for which no 2-intercept 
caustic rays emerge from a symmetric drawdown. However, an asym- 
metry of ±24 percent yields an external 2-intercept caustic ray on the 
side of reduced 7. Hereafter, the ray rotates in the forward direction, as 
do the 3-intercept rays. It is interesting to note that the / ray initially 
emerges on the reduced 7 side at ±22 percent asymmetry without 
forming an external caustic and becomes the externally visible caustic 
ray at 67 = ±24 percent. A similar ray, initially incident at the inflection 
point on the reduced 7 side as a catacaustic ray, emerges even earlier on 
the increased 7 side, at an asymmetry of 67 = ±20 percent. However, 
within the range of asymmetries considered, this ray never develops into 
an external caustic. 

Reducing the 7 to 0.688 creates a far more complicated structure even 
for the symmetric case, for which there are three 2-intercept caustics 
associated with multiple folds in the fan of emerging light. As shown 
earlier in Figs. 2a and 2c and further in Fig. 7c, the rotation is in general 
forward, although the second 2-intercept caustic appears to be relatively 
less responsive. On the decreased 7 side, the fold in the fan of emerging 
light disappears by 67 = ±10 percent, leaving the / ray still emergent. 
No first or third 2-intercept caustics are formed. At 67 = ±25 percent, 
even the / ray has disappeared. At 67 = ±30 percent, the second 
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2-intercept caustic is extinguished on the reduced 7 side and a fourth 
2-intercept appears on the opposite side. A better understanding of this 
new caustic, which is the only one that rotates backward as shown by the 
last two diagrams of Fig. 7c, may be obtained from consideration of the 
left and right diagrams of Fig. 4a. As mentioned above, the left diagram 
in Fig. 4a shows that the fan of light which lies between bounding rays 
L3 and L4 is folded to form the three 2-intercept caustics with the first 
and third very nearly superimposed. The right side of Fig. 4a shows the 
same situation on the increased 7 side for the asymmetric case repre- 
sented by Fig. 2c with the folds of the first and third 2-intercept caustics 
widely separated. The break in the 2-intercept fan of light shown lying 
between bounding rays L2 and L3 in the left diagram in Fig. 4a (which 
separates the initial fan L1-L2 from the folded fan L3-L4) is closed in 
the right diagram by the new or fourth 2-intercept caustic. This new 
caustic is folded downstream in the same way as the original first 
2-intercept caustic and originates from a coaxial ray initially propagating 
axially down the sample at an even greater radius than the second 
2-intercept caustic but well inside the ray which originates the 3-inter- 
cept caustic (see Fig. 2c). 

At a further reduction of scale to a 7 of 0.595, all rays rotate forward 
due to increasing asymmetry, with the first 2-intercept the fastest and 
the second 2-intercept the slowest. In this case, the former disappears 
for the increased 7 side at an asymmetry of almost ±30 percent, while 
the fold in the fan of light separating the first and third 2-intercept 
caustic rays on the opposite side is almost completely suppressed by an 
asymmetry of only ±20 percent. 

A final case, at a 7 of 0.450, behaves in much the same way, although 
since the 6 angles are even smaller, the first 2-intercept caustic disappears 
from the increased 7 side of the drawdown zone at an asymmetry of only 
±4 percent. In this case the 3-intercept caustic propagates at a lesser 
angle than the first 2-intercept caustic and disappears at around the 
same percent asymmetry and on the same side. 

3.2 Discussion 

Study of the collective response yields two general observations. First, 
it is apparent that all of the significant caustic rays rotate in a forward 
direction with increasing asymmetry: the first 2-intercept caustic being 
most responsive, the 3-intercept caustic slightly less so, and the third 
and second 2-intercept caustics even less so, in that order. Second, the 
angles between the opposing rays for each caustic appear to be fairly 
constant with changing asymmetry, except when one of the rays nears 
extinction by internal reflection. 

Figures 8a, 8b, and 8c illustrate these points. Each of these figures is 
a plot of the propagation angle 6 as a function of the percentage of 7 
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Fig. 8 — Plots of the 8 angles of the caustic trajectories for various 7s as functions of the 
asymmetry, f>y. Here the values of the angles, II-, emergent on the side where 7 is decreased 
are plotted on the left; the values of the angles, 0+, emergent on the side where 7 is increased 
are plotted on the right. These results are arranged as follows: (a) The first 2-intercept 
caustic angles, designated 8 UI , and the / ray propagation angle, 0', vs the asymmetry, f>-y, 
in percent, (b) The 3-intercept caustic angle, D , vs the asymmetry, 67, in percent, (c) The 
second and third 2-intercept caustic angles, 6 m and f/5 , vs the asymmetry, 67, in per- 
cent. 

asymmetry, dy. The values for the side where 7 is decreased, 0_, are 
plotted to the left; and the values for the side where 7 is increased, 0+, 
are plotted to the right of the zero asymmetry line. Figures 8a and 8b 
represent the response of the first 2-intercept and the 3-intercept caustics 
respectively, while the responses of both the second and third 2-intercept 
caustics are presented in Fig. 8c. The 6 response of the / ray is repre- 
sented by a broken curve where it appears alone on Fig. 8a, and broken 
curves are used in Fig. 8c to differentiate between the second and third 
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2-intercept caustics. Each curve shown is labeled to indicate the average 
7 for which the range of asymmetries is now determined. Thus, to de- 
termine the 8 angles for a first 2-intercept caustic ray on a melt zone of 
average 7 of 1.323 with an asymmetry of dy = ±30 percent, one finds a 
value of 115 degrees on the left side of the plot as the value of 8- 1 on the 
side of the melt zone where y is reduced 30 percent, and a value of 77.4 
degrees on the right of the plot as the value of 8+ l on the side of the melt 
zone where y is increased 30 percent. 

Except for those portions of the curves which are near points of ex- 
tinction, the curves* for both the first 2-intercept caustic and the 3- 
intercept caustic are similar in shape for all 7s. In fact, over the range 
of most realistic asymmetries, which means the first ±10 or ±15 percent, 
these curves are almost linear. This indicates that for these caustics the 
total included angle between opposing caustic rays, 8+ + 8-, remains 
constant while the angle of its bisector increases linearly with the in- 
creasing asymmetry. In terms of actual drawdown samples, these caustic 
rays form cones of light which remain fairly stable in shape but which 
tilt in the direction of these asymmetries and become distorted in those 
regions where they approach extinction by internal reflection. This is 
also true of the second 2-intercept caustic except that it seems to be 
relatively less sensitive to this asymmetry. 

Another way of illustrating this behavior is shown in Fig. 9. Here the 
average sum, 83 = %(0- + 8+), is a direct measure of the constancy of 
the cone half-angle as a function of asymmetry, while the average dif- 
ference, 8d = l k(8- — 8+), gives a direct measure of the tilt or rotation 
of the cone as a function of asymmetry. Figure 9a shows that the cone 
half-angle for the first 2-intercept caustic, Bg 1 , remains fairly constant 
but does decrease a few degrees at high asymmetries. This may be par- 
tially accounted for by the fact that the average angle does not remain 
constant with increasing asymmetry but decreases as shown in Fig. 10. 
Consequently, if the first 2-intercept caustic is somewhat more 
/3-dependent, it too may be expected to decrease, since in Part I 1 on 
symmetric drawdowns it was shown that decreasing /? decreases 8. On 
the other hand, Fig. 9b shows fairly constant 3-intercept caustic-ray cone 
angles, 8$, for all circumstances except those approaching extinction, 
as does Fig. 9c for the 83 angles of the second and third 2-intercept 
caustics. 

In all plots of Fig. 9, the curves for the average differences, 80, vs 



* Many more calculations were made to produce these plots than were illustrated by the 
diagrams of Fig. 7, with particular attention given to accurate determination of the points 
where each caustic disappeared, if it did so within the range of interest. These terminations 
are shown exactly where they occur on the plots of Fig. 8, with arrows marking a contin- 
uation where the computation was stopped for other reasons, such as reaching the limits 
of the maximum range of interest or reaching the limits of the computable range of the 
input data. 
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Fig. 9— Plots of the averaged sums, S = V2C- + "+). and the averaged differences, D 
= x k(0- - 0+), vs the asymmetry, 67, for various 7s. These results are arranged as shown 
in 9a, 9b, and 9c. (a) The first 2-intercept caustic angle's averaged sums and differences 
vs the asymmetry, ±^7, in percent. {Figs. 9b and 9c on following pages.) 

asymmetry, by, are nearly linear within the same limitation and indicate 
a constant slope or rate of rotation, rp. The greater slopes, r Dl of the 
average difference curves in Fig. 9a indicate the higher response of the 
first 2-intercept caustic to changes in asymmetry, which is about 0.8 
degree/±l percent. Also, their closer grouping indicates a significant 
degree of uniformity of this rate over a range of 7s. The 0£> vs 07 curves 
for the 3-intercept caustic ray angles are also nearly linear but shallower, 
with rates of rotation, r D , of about 0.6 degree/±l percent. They also show 
considerably more variation from 7 to 7. 

Figure 9c shows the same two functions {6s vs 67 and &d vs 67) for the 
second and third 2-intercept caustics. For both, the half-cone angles, 
0.s, are reasonably constant except near extinction (i.e., the second caustic 
at 7 = 0.68 degree), but the rates of rotation are smaller. These run from 
0.4 degree/±l percent and less for the third caustic to 0.17 degree/±l 
percent and slightly negative for the second caustic, and with much 
greater variation from 7 to 7. 

The results provide a means by which the caustic angles measured on 
opposite sides of a section through the melt zone may be used to estimate 
both average 7 and its asymmetry for that section. For example, since 
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Fig. 9(b) — The 3-intercept caustic angle's averaged sums and differences vs the asym- 
metry, ± l>y, in percent. 



the averaged sum of the 3-intercept caustic angles, 6$, is relatively in- 
sensitive to the asymmetry in 7 (Fig. 9b), this parameter may be used 
to estimate 7. In other words, 6$ will be very nearly the same as D for 
the symmetric case. In turn, the dependence on 7 may be derived directly 
from the 6 D vs relationship given in Fig. 6 of Part I. 1 Figure 11 shows 
a plot of 6g vs 7 adapted from this original plot of d D vs /3 using 7 = tan 
/3 and 0§ « 8 D . The figure also shows a plot of the inverse rate of rotation, 
r D X vs 7- Therefore, for a drawdown where melt zone asymmetry is the 
result of an asymmetry in 7, measurement of 0+ and 6- will provide a 
means of estimating 7 and its asymmetry as follows: 

From 0£ and 0?, compute 6% = W- ~ 0+) and Og = y 2 (d£ + 0?). 
Then use Fig. 11 to estimate a value of 7 and to find the corresponding 
value of ro 1 . The latter may then be used to compute the asymmetry, 
by, by multiplying by 0$- Consider the example of a sample whose 
3-intercept caustic angles measure 134 and 115 degrees. These yield a 
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Fig. 9(c) — The second and third 2-intercept caustic angles' averaged sums and differ- 
ences vs the asymmetry, ± ^7, in percent. 

ds = 124.5 degrees which, by Fig. 11, gives an estimated 7 of 0.930. At 
this 7, Fig. 11 also gives an r^ 1 = 1.56. Since Oft = 9.5 degrees, the 
asymmetry, &y, is then 1.56 X.9.5 = 14.8 percent. Computer-generated 
data for a 7 of 0.926 with an asymmetry of ±15 percent gives d_ = 134.1 
degrees and 0+ = 115.0 degrees, which is within an absolute error of only 
0.004 in estimating 7 and 0.2 percent in estimating the asymmetry. This 
order of accuracy is consistent throughout the region for which 6$ is 
independent of 7. 

The development of additional data for the various 2-intercept 
caustics would permit similar analyses of data taken over the appropriate 
7 or j3 ranges for which these caustics appear on both sides of the melt 
zone and manifest linear dependences on by or 50. 

IV. COMPARISON WITH AN ACTUAL DRAWDOWN SAMPLE 

Empirical profile data from Sample 1 discussed in Ref. 2 provide a 
comparison for the present technique. Figure 12 shows data taken from 
the 3 and 9 o'clock profiles superimposed on smooth curves computed 
from the averaged data used in the present report (from Sample 4). 
These synthesized profiles have been asymmetrically scaled (7 = 1.006 
± 15 percent) to match the angles of the empirical data (0+ = 48.9 
degrees (3-side) and /3_ = 44.6 degrees (9-side)). It can be seen that the 
empirical and computed values are very similar and superimposed (in 
this case) with no axial shift asymmetry (A/yR = 0). Various caustic 
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Fig. 10— Plots of the averaged sums, S = %(0- + 0+), and the averaged differences, 
fj D = %(/3_ — 0+), vs the asymmetry, 67, in percent, for various 7s. Here 0- is the incident 
angle on the side of reduced 7 and 0+ is the incident angle on the side of increased 7. 




Fig. 11— Plots of the averaged sum of the caustic angles, OQ, and the inverse rate of 
rotation (from their averaged difference curves), r5\ vs 7 for the 3-intercept caustics. 

angles were computed from these profiles and are shown in Table II along 
with the angles measured in the original experiment and the values 
calculated at that time using the empirical profile data. In some cases, 
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7= 1.066 ± 15%. A/7 R = 0.0 
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Fig. 12 — Plots of the empirical profile data from Sample 1 (Ref. 2) and the computed 
profiles obtained using the present technique for synthesizing asymmetric profiles (7 = 
1.006 ± 15 percent, S/yR = 0) with differing slopes (/3+ = 48.9 degrees, /3_ = 44.6 degrees). 
For Sample 1, Ro = 3.24 mm, and for Sample 4, Ro = 3.09 mm. 

Table II — Caustic angles from profile data 



Measured Values From 

Original Experiments 

(Ref. 2) (degrees) 



Calculated Values Using 

Empirical Profile Data 

(Sample 1) (degrees) 



Computed Values Using 
Smoothed Avg'd Data 
(Sample 4) (degrees) 



•J? 

0" 



46.1 

72.5 
130.6 
140.0 

59.3 
135.3 

13.2 
4.7 



54.5 

75.3 
133.8 
142.7 

64.9 
138.3 

10.4 
4.5 



EXT (43.2)* 

75.Q 
122.1 
145.4 
(59.1)* 
133.8 
(15.9)* 

11.7 



* In this case the caustic is extinguished by internal reflection, so the extinction angle itself 
is used in computing effective 0^ and flp values. 

the newer values are in better agreement with the observed angles than 
are those computed using the actual profiles, and certainly they appear 
to be about as good overall. 

V. SUMMARY 

Two types of asymmetry have been studied using computer simulation 
for predicting the trajectories of caustics emitted by a coaxially illumi- 
nated, fused-silica melt zone. The first type of asymmetry involves an 
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axial shift of one side of the melt zone with respect to the other. Increases 
in this type of asymmetry produce changes in the trajectories of the 
various types of caustics which appear. The appearance of these caustics 
also depends very strongly upon the average taper of the melt zone. 2 The 
following effects have been observed: 

(i) The first 2- and 3-intercept caustics rotate downstream when 

they appear on the trailing side of the melt zone. 
(ii) The first 2-intercept caustic initially starts to rotate upstream 
and then reverses direction when it appears on the leading side. 
This effect is more pronounced in blunter melt zones. 
(Hi) The first 3-intercept caustic generally rotates upstream on the 
leading side but may reverse its rotation for large asymmetries 
on very blunt melt zones. 
(iu) The second 3-intercept caustic appears on the leading side only 
for blunt melt zones (7 = 0.450) and only at large asymmetries 
for which it rotates downstream. 
(v) The second and third 2-intercept caustics rotate downstream 
on the leading side and upstream on the trailing side with in- 
creasing asymmetry. This is the reverse rotation of the first two 
major caustics. 
(ui) The fourth 2-intercept caustic appears only on the leading side 
of melt zones which are moderately blunt (7 < 0.688) and ro- 
tates downstream with increasing asymmetry. 
{vii) When a particular caustic appears on opposite sides of a melt 
zone, its intersecting trajectories form an angle which may be 
thought of as the included angle of a cone of light emitted from a 
three-dimensional sample.* Cone angles of the first 3-inter- 
cept caustic and the second and third 2-intercept caustics are 
moderately insensitive to the asymmetry. 
(uiii) The cone angle of the first 2-intercept caustic increases signif- 
icantly with increasing asymmetry. 
{ix) The cones of light for most of these caustics rotate forward with 
increasing asymmetry. The rotation of the cones of light for the 
3-intercept caustic and the first 2-intercept caustic is initially 
forward or toward the leading side and then reverses. 
(x) The light cones for the second and third 2-intercept caustics 
generally rotate backward at a constant rate with increasing 
asymmetry for a given 7. 
(xi) The rate of rotation of the first 3-intercept caustic is generally 
greater than that of the first 2-intercept caustic, while the rates 
of rotation for the second and third 2-intercept caustics are al- 
ways negative. 



* It should be kept in mind that the cone angle, as we have defined it here, is the sum of 
the particular caustic angles on opposite sides of the drawdown measured from the up- 
stream axial direction. 
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The second type of asymmetry investigated involves scaling opposite 
sides of the melt zone differently to produce asymmetries in slope. In- 
creases in slope asymmetry produce changes in the external caustic 
trajectories, which are in many ways simpler than those observed with 
the axial shift asymmetry. These changes are as follows. 

(0 All caustics generally rotate downstream on the trailing (de- 
creased) 7 side and upstream on the leading (increased) 7 side. 
This is comparable to the behavior of the first 2- and 3-intercept 
caustic response to the shift asymmetry. 
(ii) The cone angles of the 3-intercept caustic and the second and 
third 2-intercept caustics are almost independent of the asym- 
metry over significant ranges. 
(Hi) The cone angle of the first 2-intercept caustic is only weakly 
dependent on the asymmetry in 7 and may be varying in response 
to changes in the average value of /3, which itself changes with 
asymmetry. 
(iv) There is no second 3-intercept caustic, as was the case previ- 
ously. 
(u) All the caustics rotate forward at constant rates for any partic- 
ular 7. 
(ui) The rates of rotation are ordered as follows: 

I'm > r D > r U3 > r U2- 

For both types of asymmetry, sufficient rotation in any axial plane 
usually causes one of the caustic trajectories to be extinguished by in- 
ternal reflection. The approach of a caustic trajectory to its extinction 
usually results in an acceleration in its rate of rotation, thereby distorting 
the general pattern of movement of the cone of light for that caustic. 
Nevertheless, over much of the range of reasonable geometries there will 
be caustics whose rotations are sufficiently regular as to provide a direct 
means of estimating the geometrical asymmetry from the asymmetry 
of the caustic cone itself, if the caustic and type of asymmetry are 
properly identified. 
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